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■ Abstract. We examine the binding conditions for atoms in non-relativistic QED, and 

prove that removing one electron from an atom requires a positive energy. As an application, 
we establish the existence of a ground state for the Helium atom. 
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1. Introduction 

One of the most fundamental results in the spectral theory of multiparticle Schrodinger 
operators is the proof of the existence of a ground state for atoms and positive ions. It was 
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accomplished for the Helium atom by T. Kato in 1951 jU], and for an arbitrary atom by 
G. Zhislin in 1960 [12J (cf. the Zhislin theorem in [TT]). 



The standard approach to the proof of these results consists of two main parts. The 
first key ingredient is the HVZ - (Hunziker - van- Winter - Zhislin) theorem, which establishes 
the location of the essential spectrum, and gives a variational criterion for the existence of 



a bound state. The latter can be referred to as "binding conditions". The statement is that 
the bottom of the essential spectrum of the whole system is defined by its decomposition 
into two clusters. If the infimum of the spectrum of the entire system is, for all nontrivial 
cluster decompositions, less than the sum of the infima of the spectra of the subsystems, it 
follows that the whole system possesses a ground state. 

For an atom with infinite nuclear mass, this condition can be written as 



(1) E V {N) < E V (N') + E°{N - N') for all N' < N, 

where E V (N) is the infimum of the spectrum of the atom, E V (N') is the infimum of the 
spectrum of the same atom without (N — N') electrons, and E°(N — N') is the infimum of 
the spectrum of the system of (N — N') electrons, which do not interact with the nucleus. 
Obviously, in the case of Schrodinger operators (in Quantum mechanics) E°(N — N') = 0, 
and according to the HVZ theorem, it suffices to consider only the decompositions with 
N' = N- 1 in HJ). 

The second key ingredient consists of the construction of a trial state for the Hamil- 
tonian of the whole atom with energy less than E V (N — 1). As noted above, this step was 
accomplished by T. Kato for Helium, and by G. Zhislin for the general case. 

The problem of the existence of the ground states of atoms has attracted new atten- 
tion in the context of non-relativistic quantum electrodynamics in the more recent literature. 
Bach, Frohlich and Sigal [2] first established the existence of the ground state for the ultra- 
violet regularized Pauli-Fierz Hamiltonian of an atom, for sufficiently small values of some 
constants in the theory. 
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It was subsequently established in [S] that the criterion for the existence of the ground 
state of multiparticle Schrodinger operators can be extended to hold for Pauli-Fierz Hamil- 
tonians in non-relativistic QED, for arbitrary values of the parameters of the theory. 1 

The problem, however, of devising a mathematically rigorous proof of the fact that 
the binding conditions are fulfilled for atoms apart from the one electron case, which was 
covered by jHj, has turned out to be very complicated. To clarify the main obstacles, let us 
recall the basic idea underlying the proofs of the Kato and Zhislin theorems. 

If the system is separated into a pair of clusters, one of which contains N —1 electrons 
close to the nucleus, and the other comprises a single electron far away, there is an attractive 
Coulomb potential that acts on the separated particle. If the latter is localized in a ball 
of radius R centered at some point with distance bR from the origin, and the subsystem 
with N — 1 electrons is localized in a ball of radius R centered at the origin, the intercluster 
Coulomb interaction can be estimated as CR~ l with C < for b > N. At the same time, 
localizing the subsystems in these balls requires an energy CR~ 2 in the case of Schrodinger 
operators. For large R, the Coulomb term is obviously dominant, and the binding condition 
is fulfilled. 

This is contrasted by the situation in non-relativistic QED, where the particles have 
to be localized together with the quantized radiation field. One can expect, on the basis of 
dimensional analysis [S], that such a localization requires an energy CR~ l , which makes it 
impossible to establish the dominance of the Coulomb interaction by scaling arguments. 

In the work at hand, it is demonstrated how this obstacle can be overcome. We prove 
that if the self-energy operator T , restricted to states with total momentum 0, possesses a 
ground state, it is possible to construct a state consisting of an electron coupled to a photon 
field, localized in a ball of radius R with energy Eo + o(_R -1 ), where So is the self-energy of 
an electron. Hence, similarly as for Schrodinger operators, the localization term o(i? _1 ) can 
again be compensated by the attractive Coulomb potential. This implies that the binding 
condition is fulfilled for decompositions into clusters with N — 1 and 1 particles. 

Existence of the ground state of T has been recently established for sufficiently small 
values of the fine structure constant |3]. It was proved earlier in [H] that for the decomposition 
into clusters with zero electrons and N electrons, the binding condition is also fulfilled. Thus, 
if an atom or a positive ion has only two electrons, the ground state exists. 

If an atom has more than two electrons, one must also verify the binding conditions 
for 1 < iV — N' < N. We note that in contrast to the quantum mechanical case, a system 
of K electrons coupled to a photon field may have an energy smaller than the self-energy of 
an electron multiplied by a factor K. 

To control this case, it would be sufficient to combine a straightforward modification of 
the method developed in this paper with a generalization of the results of [S], and to apply 
it to the case of a system without external potential, after separating the center of mass 
motion. This generalization is, however, beyond the scope of the present work. 

The first proof of the existence of the ground states for all atoms in non-relativistic 
QED has, besides numerous other important results, been accomplished by Bach, Frohlich 
and Sigal in , by a completely different approach. To compare the results in j^j for Helium 
to the results of the work at hand, we remark that the units used in our paper correspond 
to those in jH], which differ from the ones in [2]. Furthermore, we emphasize that while 



A detailed review of numerous further results connected to the existence of ground states, mostly in 
Nelson- type models, can be found in [S]- Furthermore, also cf. [7| 
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the ultraviolet cutoff in the quantized vector potential employed in [2] is, in our units, 
incorporated at a value A ~ a, where a denotes the fine structure constant, we are studying 
the corresponding case for an ultraviolet cutoff at A ~ 1. The parameter that accounts for 
the strength of the perturbation produced by the photon field is in assumed to be much 
smaller than a constant that depends on the ionization energy of the atom, the latter being 
computed for the Schrodinger operator of the electron subsystem. One of the key issues in 
the work at hand is to devise a proof that also encompasses the strongly nonperturbative 
regime, where this parameter is allowed to be much larger than the ionization energy. This 
is achieved mainly based on the parameter independence of the results of jH], as well as of 
the methods developed in the present paper, in addition to exploiting the existence of the 
ground state of T for small a. 

2. Definitions and main results 

We consider the Pauli-Fierz Hamiltonian Hn for a system of A" electrons in an external 
electrostatic potential, coupled to the quantized electromagnetic radiation field, 

JV 



The operator Hn acts on the Hilbert space H := H N ® T, where H N , for JV < oo, is 
the Hilbert space of A" non-relativistic electrons, given by the totally antisymmetric wave 
functions in (L 2 (IR 3 ) ® C 2 )^, where IR 3 is the configuration space of a single electron, and 
C 2 accomodates its spin. 

We will describe the quantized electromagnetic field by use of the Coulomb gauge 
condition. Accordingly, the one-photon Hilbert space is given by L 2 (M 3 ) <g> C 2 , where M 3 
denotes either the photon momentum or configuration space, and C 2 accounts for the two 
independent transversal polarizations of the photon. The photon Fock space is then defined 



where the n-photons space Ti = (-^ 2 (1R 3 ) ® C 2 ) is the symmetric tensor product of n 
copies of L 2 (M 3 ) <g>C 2 . 

We use units such that H — c—1, and where the mass of the electron equals m — 1/2. 
The electron charge is then given by e = y/a, with a ~ 1/137 denoting the fine structure 
constant. As usual, we will consider a as a parameter. 

The operator that couples an electron to the quantized vector potential is given by 



where by the Coulomb gauge condition, divAf = 0. The operators a x , a* x satisfy the usual 
commutation relations 



(2) 




l<k,£<N 



by 




e x {k) e ikx ® a x {k) + e~ ikx ® a\{k) Ak =: D{x) + D*{x), 



[a u (k),a* x (k')]=5(k-k')6 X! 



[a v {k),a x {k')} = 0, 
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and there exists a unique unit ray Qf G J 7 , the Fock vacuum, which satisfies a\(k)Qf = 
for all fceR 3 and A € {1,2}. The vectors £\(k) € M 3 are the two orthonormal polarization 
vectors perpendicular to k, 

= —===- and £ 2 {k) = — Ae^k). 
y/kf + k\ \k\ 

The function C(l^l) describes the ultraviolet cutoff on the wavenumbers k. We assume ( to 
be of class C 1 , with compact support. 

The operator that couples an electron to the magnetic field Bf = curl A f is given by 



X=l,2 J 1 1 

=: K{x) + K*(x). 

In Equation ([%]). a = (eri, o~2, 03) is the 3-component vector of Pauli matrices 

o i\ /o -A A o 



dA; 



ai= »1 ' U ' ^ 10 -1 



The photon field energy operator Hf is given by 

H f = J2 [ \k\a* x (k)a x (k)dk. 

X=1,2 J ^ 

The potentials V and W are relatively —A bounded with relative bound zero and 
satisfy for positive 7, 70 and r the following conditions: 

(3) V{x) < \x\ > r , 

\x\ 

(4) W(x) < ^r, \x\ > r . 

One of the main assumptions of the work at hand is the existence of a ground state of the 
one electron self-energy operator with total momentum P = 0. For its precise formulation, 
let us consider the case of a free electron coupled to the quantized electromagnetic field. The 
self-energy operator T is given by 

T = (-iV x ® If + v 7 "^/^)) 2 + v 7 "^ • + ^ ® 

We note that this system is translationally invariant, that is, T commutes with the operator 
of total momentum 

Plot = Pel® If + hi® Pf, 

where p e i and Pf = J2x=i 2 I ka* x (k)a\(k)dk denote the electron and the photon momentum 
operators. 

Let Hp = C 2 ® T denotes the fibre Hilbert space corresponding to conserved total 
momentum P. For any fixed value P of the total momentum, the restriction of T to the 
fibre space Hp is given by (see e.g. [3]) 

(5) T(P) = (P-P f + ^A f (0)) 2 + y^a ■ 73/(0) + H f . 

We denote £ = inf a(T) and E = inf a(T(0)). The following assumptions will be used to 
formulate the main result 
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Condition £ L . 

i) S = S 

ii) S is an eigenvalue of T(0), with associated eigenspace £s . 

iii) There exists Qq £ £e with a finite expectation number of photons, i.e. 

(N f n ,n ) < c, 

where Nf = J2 x =i 2 I a*\(k)a\(k)dk. 

iv) The above eigenfunction Q fulfills, for A = 1, 2 and some po G (6/5, 2] 

||V fc a A (A;)fio|| G^ P0 (K 3 ) + L 2 (M 3 ) 

Condition i) was studied by Frohlich for a spinless Pauli-Fierz model, jB], who proved that 
in this case, it is fulfilled for all a > 0. 

For the case including the a -B term, it was proved in ^ that for small a, the condition 
is also fulfilled. 

The existence of the eigenspace S-^ in ii) was recently proved for sufficiently small a 

Finally, it will be proved in the present paper that for small a, the function Qq possesses 
the properties iii) and iv). Thus, we conclude that there exists a number ao, such that at 
least for all a < ao, condition <£ L is fulfilled. 

The second main set of assumptions required for our analysis is given as follows. For 
M G N, let Hm denote the Pauli-Fierz Hamiltonian for M electrons defined in (0). 

Condition (£ 2 - 

i) The operator Hm has a ground state 

(6) t e n = n% ® t, 

with a finite expectation number of photons. 

ii) For A = 1, 2 and some p G (6/5, 2], 

||(/ei® V fc a A (fc))T|| G L po (M 3 ) +L 2 (M 3 ). 

iii) Let Xi i = 1, . . . M be the position vectors of the electrons. Then, 

M \ 

^2\xi\®I f TeW. 
i=i / 

For M G N, let 

E M = Ma(H M )- 
The main result of this article is the following 

Theorem 2.1. For N G N, /et i/ie Conditions C ± and C 2 M = N — 1 be fulfilled, and 
assume that the potentials V and W satisfy (jSJ) and (JU), with 70/71 > (N — 1). T/ien, 

(7) .Etv < £jv-i + S. 

Remark 2.1. If one assumes that the system with M electrons satisfies the binding condition 
°f 0; it was shown in [8 j t/iat this system possesses a ground state which satisfies all the 
conditions 0/C2. In particular, the ground state of the Hydrogen atom fulfills £2- 
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This Theorem shows that under the above stated conditions, removing one electron 
from the system costs energy. In this sense, the system is stable with respect to the given 
type of ionization. 

The conditions on the potential V(x) and W(x) cover a large number of models in 
atomic and molecular physics. In particular, for V(x) = —(3Z/\x\ and W = /3/\x\, the 
operator Hn describes an atom or ion with N electrons. 

In the physical case, f3 is equal to the Sommerfeld fine structure constant a. However, 
we would like to emphasize that the proof of the Theorem is valid for all values of (3 > 0, 
even in the strongly nonperturbative regime < ,3 < a. 

Theorem 12.11 states that as long as the number of electrons N is less than Z + 1 
(neutral atoms and positive ions), ionization by separation of one electron is energetically 
disadvant ageous . 

If was earlier proved in jH] that removal of all electrons from the atom also leads to an 
increase of the energy. 

Combining these two results for the case N = 2, and the binding condition in jHJ 
Theorem 3.1], yields 

Theorem 2.2. The Pauli-Fierz Hamiltonian for Helium 

2 

H 2 =J2 {("^ ® J / + V^Mx*)? + >far ■ B f (xt) - g ® I/} 
i=\ 

has a ground state for all a < ciq. 

Notice that the conditions on the potential V(x) require only some type of behaviour 
at infinity. Therefore, instead of one nucleus with Coulomb potential of charge Z, one can 
consider a system of nuclei 

V(x) = £ j^- 

i=l ' l ' 

with the same total charge, in the infinite mass approximation. In particular, for Hydrogen 
molecules as well as for all molecular ions with two electrons, Theorem 12.11 implies the 
existence of a ground state for all a < ao- 



3. Properties of the ground state of T(0). 

This section addresses the main properties of the self-energy operator T(0) that are 
required for the present analysis. In particular, existence of a ground state fl G C 2 <E> J 7 , 
finiteness of the expected photon number with respect to Q , and regularity of ax(k)Q are 
discussed. 

3.1. Existence Theorem. In the following theorem, existence of a ground state of T(0), 
and bounds on the associated expected photon kinetic energy are established. 

Theorem 3.1. For a sufficiently small, So = infcr(T(0)) is a degenerate eigenvalue, border- 
ing to absolutely continuous spectrum, which satisfies 

| S 1 < ca . 
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Let £s = ker(T(0) — S ) C C 2 <8> T denote its eigenspace. Then, dim c £ So = 2, and for any 
£ £e ? normalized by (f2 ,fi/) = 1, the estimate 

||^o || < 1 + c\/a 

satisfied. Furthermore, 

(8) P/(0)^o|| , ||^) /2 ^o|| <cv^ 

/jo/(i. t4// constants are uniform in a. 

For the spinless case, both results are proved in jH] by use of the operator-theoretic 
renormalization group based on the smooth Feshbach map, cf. pQ. For the case including 
spin, an outline of the proof is given in the Appendix of 4J, while a publication containing 
the detailed proof is in preparation. The bound on 1 1 ^4 y (0) 1 1 follows straightforwardly from 
the one on ||ifyf2 ||. 

3.1.1. Expected photon number. Using Theorem l3.ll we may next bound the expected photon 
number with respect to no- 



Theorem 3.2. For a sufficiently small, and Qq G £e defined as in Theorem IJ.il Qq G 
Dom(A^y 2 ) ; where Nf = Ylx=i 2 I ct>\(k)a x (k)dk is the photon number operator, and 

HiV^oll 2 < cyfa. 

In particular, 

\\xi\k\<\)a x ik)^\\<c^\k\- x . 
All constants are uniform in a. 

Proof. We first remark that the integral J dk\\a x (k)Qo\\ 2 is ultraviolet finite, since 
X (\k\ > l)||a A (A;)fio|| 2 dA; < J X {\k\ > l)\k\ \\a x (k)Q \\ 2 dk 

(9) < ca , 

using (|HJ). We may thus assume that the domain of the integral is the unit ball -Bi(O). For 
I A; I < 1, we employ a similar argument as in [HIEllHj. Using 

( : T(0) : -Z' )a x (k){l = [ : T(0) : , a x (k)]Q , 

where : ( • ) : denotes Wick ordering, and 

E' :=E -(A f (Of) Qf = mia(:T(0):), 

we obtain 

a x (k)n = ^R(h)(k-A f {0) + ^^e x (k)-P f 

(10) + TkW ik A ex{k) ■ ° + ^Jkjjl 6 ^ • A /(°)) fi o > 

where 

(11) R(k) := (H, + \k\ + \{P f + kf - S' ) _1 . 



Clearly, (Qf, : T(0) : Qf) = 0, and a standard variational argument shows that £(, < for 
a > 0. Hence, < R(k) < (H f + \k\y\ and 

\\R(k)P,\\ < \\R(k)H f \\ < 1 . 

Thus, using < 1 and theorem 13.11 

||x(|fc| < l)a A (A0fto|| < c v ^ X (|A:|<l)(p / (0)fio||+2|A;|- 1 / 2 ||fio|| 

+ v / ^|A;r 1 p / (0)fio| 

(12) < Cy /a\k\- x . 

The right hand side is in L 2 (B>i(0)), and the assertion is established. □ 

For the case of a confined electron, it was proved in jH] that the corresponding estimate 
exhibits a l^l" 1 / 2 singularity instead of as present here, owing to the exponential decay 
of the particle wave function. 

Furthermore, if the conserved momentum P is non-zero, there exists a ground state 
Qp(k) for a regularized version of the model, which includes an infrared cutoff below < 
k 1 in v4j(0) (some requirements on the cutoff function are necessary, cf. jS])- Then, with 
all modifications implemented, the additional term 

^R(k)^$P ■ e x (k) n P (K) 

enters the right hand side of (jlOj) . Therefore, (flp(/c), NfQp(n)) is logarithmically infrared 
divergent in the limit k — > 0, for all \P\ > 0, and in fact, Qp(k) does not converge to an 
element in Fock space. 

3.1.2. Regularity properties of the ground state. Next, we derive a result about the regu- 
larity of a\(k)Q in momentum space, which is, in our further discussion, used for photon 
localization estimates in position space. 

Theorem 3.3. For a sufficiently small, let Q e £s - Then, 

||V fe a A (A;)no|| GL p (M 3 ) + L 2 (M 3 ) , 

for 1 < p < | . 

Proof. We proceed similarly as in jSj. To begin with, we differentiate the right hand side of 
(fTUj) with respect to k, and observe that 

(13) \V k R(k)\ < (l + H f + \k\)R 2 (k) , 
since \Pf\ < Hf. 

Let us first bound the ultraviolet part of || V^a-A^^oll- For \k\ > 1, 
||x(|A:| > l)V*flA(A;)no|| = v^llxd^l > l)V k R{k)k ■ A f {0)tt \\ 

< (\\ x (\k\>l)(l + H f + \k\)R(k)\\ + \k\ 
\\V^ X (\k\>l)R(k)k-A f (0)Q \\ 

(14) < 2v^||x(|A:| > l)ax(k)n \\ , 



and consequently, by Theorem 13.21 

(15) / \\V k a x (k)Q \\ 2 dk < ca . 

J\k\>i 

We may thus restrict our discussion to the case \k\ < 1. 

Differentiating with respect to k, the photon polarization vectors satisfy 

c 



(16) |V - WI £ MTW 

Recalling that the cutoff function ( is of class C 1 , and using Theorem 13.21 one straightfor- 
wardly deduces that there exists a constant c which is uniform in a, such that 

1 1 



< l)V fc a A (A;)fio|| < + 



(17) < ^ 



\k\^/W+M' 

Here, one again uses ||i?(A;)P/|| < \\R(k)Hf\\ < 1, and ||i?(A;)|A;||| < 1, in addition to (H3J). 
Thus, by the Holder inequality, 

\\v k a x (k)n rdk" 1,r 



\k\<l 



1/r* 

dk' 



'\k\<l \ K \ ' l K l + K 2> ' v ^|fe|<l \ K \ 

with - = i + pr. The integrals on the right hand side of (fTHj) are bounded for the choices 

1 < r* < 6, and 1 < r < 2, which implies that 1 < p < |, corresponding to the exponent 
expected from scaling. □ 

In the case of a confined electron, IHl, the bound analogous to (O is — , cv ^ „. The 

reason for the fact that it is by a factor Ik] 1 / 2 less singular is stated in a previous remark. 
Consequently, in |HJ, the inequality corresponding to (HHJ) likewise requires the choice r < 2, 
but in contrast, r* can be chosen arbitrarily large. Therefore, the result proved in [8J holds 
for i > i + J- = i that is, 1 < p < 2. 

4. Self-energy of localized states with total momentum P = 

The goal of this chapter is to arrive at a sharp upper bound on the infimum of the 
quadratic form of the operator T(0), when restricted to states where all photons are localized 
in a ball of radius R centered at the origin. 

To this end, we recall that for the Schrodinger operator —A corresponding to a free 
electron, the infimum of the spectrum on the whole space is zero, whereas the infimum on 
functions supported in a ball of radius R, with Dirichlet boundary conditions, is C/R 2 . 

The main result of this section is the following. 

Theorem 4.1. For all R > 0, there exists a function $ fl G 2D(T(0)) ; such that 
i) The n photonic components $^(yi, . . . ,y n ; X±, . . . , X n ) fulfill 

supp$^ C {{yi, . . . , y n ; A 1; . . . , A n ) | sup \y t \ < R} 

i 
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ii) 

(19) (T(0)$ K , < (s + ||$ H || 2 
where c(R) tends to zero as R tends to infinity. 

iii)The function $ R has the following additional properties. For all e > and all \x\ > 2R, 

(20) \(D(x)$ R , $ R )\ < ^j^-\m 2 , 

(21) K^V^K^l^lpi^ 

(22) < 4Sll $ l| 2 > 



\x\ 



and 



(23) UK(x)$ R ,$ R )\<?^\\M 2 

\x\ 

where c(\x\) tends to zero, uniformly in R, as \x\ tends to infinity. 

Before addressing the proof of Theorem 14.11 we shall first demonstrate how it can be 
employed to construct a state in Hi ® T that accounts for a system consisting of an electron 
coupled to a photonic field, localized in a ball of radius R centered at a fixed point b, with 
energy close to the self-energy S . For that purpose, let us, for given x G M 3 , define the shift 
operator r x : T — > T , which, for = (0 O , 1; . . . , n , . . .) G J 7 , is given by 

T x <f> n (yi, ■ ■ ■ , Vn, Ai, • • • , A n ) = 4> n (yi -x,...,y n -x;X 1 ,...,X n ). 

Theorem 4.2. Lei f be a real valued function in Cg(lR 3 ) (g> C 2 , supported in the unit ball 
centered at the origin. For R > and b G R 3 , we define Q Rb G 7Yi ® F by 

Then, for all e > and i? Zarge enough independent of b, we have 

(25) (((iV* ® // + v^/Oz)) 2 + V^a.Bf(x) + I e i ® H f ) Q R \ Q Rb ) < S + ^. 



Proof of Theorem For a real valued function /, let f R,b (x) := f(x/R — b). Obviously, 



(26) 

According to Theorem 14.11 the second term on the right hand side can be estimated by 

\\f\\ 2 (T(0)****) cCR) 

\\f R > b (x) <g> $^|| 2 - R ' 
10 



For the first term on the right hand side of (|26|) . we have 

\\f R ' b (xW ~R 2) 

which completes the proof of the Theorem. 

4.1. Localization estimates. In order to prove Theorem 14. 1[ we consider the ground state 
Qo of the self-energy operator T(0) at zero momentum, and act on it with two spatial 
localization functions U n and V n , which constitute a partition of unity {U n ) 2 + (V n ) 2 = 1 
on T. This yields a state for which all photons are inside the ball of radius R, and another 
state for which all photons are outside the ball of radius R/2. 

Clearly, the expectation of T(0) with respect to Qq is not equal to the sum of the 
expectations with respect to the two localized states. The difference, which is usually called 
the localization error, must be estimated to obtain an upper bound on the self-energy of the 
localized state. In the present subsection, we estimate the localization errors for different 
terms in the operator T(0). 

Let us to begin with define spatial cutoff functions u and v as follows. We pick u £ 
C °°(M + ) such that 



(27) u(x) 



1 if x £ [0, 1/2] 
if x > 1 



< u < 1 and v := y/1 - u 2 £ C 2 (R + ). 

For Y = (y%,y2, ■ ■ ■ ,Un) £ K. n j w e denote H^Hoo = niax!<j< n \yi\. For n £ N and all 
Y £ lR n , we also define u*(Y) = u( J! % £ ) and v^(Y) = y/1 - u*(Y) 2 . 

Next, we introduce a pair of operators lA n and V n on JF by 

(28) U R ^ = (^0, Uf (yi)MVl), • • • > «n ((Vh ■ ■ ■ VM{V1, ■■■Vn):--) 

and 

(29) = (Vo, (yi)Mvi), • • • > <iiVu ■ ■ ■ vMiyi, ...y n ),---), 
where we have omitted the polarization indices from the notation. 

4.1.1. Localization error for the field energy Hf. 

Lemma 4.1. There exists c < oo such that for all e > 0, and all R large enough, 

(30) (H f U R ^,U R ^) + (H f V R ^,V R ij)-(H f ij,ij)<(N^,ij)(^+ ' W^' 2 ^ 1 



R eR 
holds for i) £ Q(H f ) n Q(N f ) . 

Proof. Since Hf maps each n-photon sector of the Fock space T into itself, it suffices to 
estimate the localization error for the n-photon component of ip. Furthermore, since Hf acts 

on a function in as nlV^J, the statement of the Lemma follows straightforwardly from 
Lemma Ol □ 

Lemma 4.2. There exists c < oo such that for all e > 0, all R large enough, 

(MuM^uMw + <|VK^)0X^r)0> - (|V|0,0) 

/ Q1 X it it it it 

[ j <(-+ c llMlyl >R)r 



R eR 
li 



holds for all(f)eC™(R 3 ). 

Proof. By [TUJ Theorem 9], we have 



(32) 



Let us consider 
(33) 



2n 2 J |y-z| 4 

|2/-^| 4 



R 



R 



Av\\ M\ 



dydz. 



dydz. 



The term with the function v can be estimated similarly. By symmetry, it suffices to estimate 
this integral in the region where \y\ < \z\. We split the integral / into three parts I\, I2, and 
I 3 , respectively, corresponding to the regions 1Zi = {\z\ < R/2}, 1Z 2 = {\z\ > R/2, \y — z\ > 
R/4} and TZ 3 = {\z\ > R/2, \y - z\ < R/A}. 

Since \y\ < \z\, we have, in the region IZi, \y\ < \z\ < R/2. Thus, in 1Z±, we have 
u{ l 4) - u{ l -f) = 0. Therefore, 



Now, for all e > 



(34) 



h < e 



R-2 



\y - A 



-dydz + 



<c[e 



R 



+ 



eR 



I0(^)| 2 

£ Jn 2 \y-A 4 

z\>R/2)f 



dydz 



where c is a constant independent of e. 

Finally, since the derivative of u is bounded, we have the inequality \u{\y\/R) 
u{\z\/R)\ 2 < c\y - z\ 2 /R 2 . This implies 



h < c 



\<Ky)\\<K*)\\v-z\ 



n 3 



\y- A 



R 2 



-dydz 



(35) 



< 



R 2 



e\<Kv)\ 2 + (l/e)\<Kz)\< 



■R-3 



\y - A 



dydz 



CP 1 

<-U\\ 2 + —Uxi\A>Rm\\ 2 - 



□ 



4.1.2. Localization error for the operator P 2 . 

Lemma 4.3. There exists c < 00 such that for all e > and all R large enough, 

(36) (p]u r ^m r ^) + (P]v R ^y R ^) - (P}M) < ^(N f M) 

holds for ip G £2(H f ) n Q(N f ) . 
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Proof. The operator Pj maps each n-photon sector into itself. Therefore, it is sufficient to 
restrict the proof to J-\ . We have 



(37) V 

+ 2^((ViO(V^n),^ n ) + ((V^)(V^n),^n) 

i,j 

Since (w^) 2 + (v^) 2 — 1, the first term on the right hand side of ()37j) is zero. Similarly, by 
rewriting the second term as 

E« V ^) 2 )(V^n),^> + ((V^XV^n), V>n) 

we find that it is also zero. Next, we note that VjVjW^ = and V;Vji^ = if i ^ j, 
because the functions u and v depend only on the ||.||oo norm. 
Thus, we obtain 

i,j 

(38) i 

where in the last inequality, we used that for some constant c, we have |Am^| < cR~ 2 and 
|Au*| < cR- 2 . □ 

4.1.3. Localization error for PfAf(0) . 

Lemma 4.4. Letij) G £}(P/A/(0)) flS)(P/) nS(iV/) ; and assume that for some p G (6/5,2], 

||V fc a A (A;)^||^GL^(R 3 ) + L 2 (M 2 ). 

Then, the inequality 

(39) 1(^/(0^,^) + (P/A/(0)V^,V^) - (P f A f (0)1>,1>) 
holds with 5 = (p — 6/5)/2. 

Proof. Throughout this proof, we will write J dy for integration over the y variable, and 
summation over the polarization A. Here and in the rest of the paper, we define G\(x) as 
the Fourier transform of the vector function 



c 

< 



R^ 5 



I A' I 



2 
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In addition, everywhere where it does not lead to any misunderstanding, we will omit the 
photon polarization index A. 
We have 



(40) 



{P f D(0)U H ^M tp) + (P f D(0)V H ip,V H ip) - {P f D(0)iP,iP) 

n [~ i=l 

+ J G(-y n+1 )iJ n+1 ~i^ U n+l^i U n + V n+l^i V n^j d Vl ■ • • d Vn+l j 

=: J^(a„ + b n ) . 



or 

,-R 



We first estimate the term a n . We denote F = u^ +1 u^ + v^ +1 v^ — l. For |y n +i| < -R/2, either 
\y n+ i\ and then u£ +1 Q0 = u^(y 1 , ...,?/„) = 1 and t^ +1 (F) = v*(y l7 ...,y n ) = 1, 
= \y k \, for some fc ^ n + 1, and then u£ fl (y) = u^i, . . .,y n ) and u^ +1 (Y) = 
v„(yi, . . . , y n ). In both cases, we get F = 0. Thus for 5 > sufficiently small, we have 

» n 

Vn + 1 G(-y n+1 )ip n+i y2 (^MFdyi ■ ■ ■ tyn+i 

J\y n+1 \>R/2 i=1 



< Vn + 1 



(41) 



(1 + \y n+1 \) l - s \G(-y n+1 )\^ n+1 \(l + \y n+1 \) 



2,S 



\y n+1 \>R/2 

T ^(P f ip) n dy 1 . . .dy n+ i 



<- 



v^TTl^n+iKi + \y n+ i\) 25 (i + \y n +i\) 



1-6 



\G(-y n+1 )\\(Pfil)) n \dy 1 ...dy n+1 
Applying the Schwarz inequality, we arrive at 



il+<5 



(42) 



<^rallv^+l^ + i(i + bn + i|) 2 1^ +1 11(1 + \y n+ i\y- 6 G\\ LHdyn+l) 

\\{P^)n\\Ll, 



where for brevity, L\ := L 2 {dy\, . . . ,dy k ). According to Lemma 1711 in the Appendix, one 
finds that ||(1 + |y n +i|) 1_5 GlU 2 (ds/ n+ i) is finite. Therefore, 

(43) E W - (llvVTT^„ + i(l + \y n+1 \) 25 \\ 2 Ll+i + \\{P^) n \\li 



We note that 



\\V k a x {k)^\\j: G L P0 (M 3 ,d£;) + L 2 (R s ,dk) 



implies 



J> + l)IIVWi (y, + M) 2 G ^ 0/2 (M 3 , dy) + L 1 ^ 3 , dy) , 
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with ^ + ^ = 1, by the Hausdorff- Young inequality. Consequently, one can straightforwardly 
verify that for S = (p — 6/5)/2, 

(44) T l|V^ + TVn+l(l + \yn+l\ni\, < C • 

f * n + l 

n 

Moreover, 

(45) ElK P /^fe <c > 

n 

since ^ G 3)(P/). Inequalities imply that 

( 46 ) EW^F?" 



Let us turn to the estimate of b n . If maxj = i r .. !n \yi\ ^ |y n +i|, then 

n 

E( M n+i(l/i> • • .2/n+i)ViM* (yi, ...?/„)+ v^iyt, . ..y n+l )ViV*(y u . ..y n )j =0 . 
i=i 

If max{yi, . . . , y n+ i} = |y n +i|) then VjW^ = Vjf^ = for all . . . y n ), such that one finds 
maxfc=i,..,n 1 I > This means that except on a set of measure zero in IR n , the functions 
u n+i^ 'i u n + V n+i^i v n have disjoint supports. Therefore, 

n 

i=i 

Moreover, V«m^ and ViV§ have support in the set {\yi\ G [R/2,R]}, thus, since from the 
above, we only have to consider the region where |y n +i| > m a x i=i,...,n \yi\, we get |y n +i| > -R/2, 
hence 

\b n \ < ^Vn / |G(-y n+ i)| |V>n+i| \i>n\dyi ■ ■ • dy«+i 

X J\y n+1 \>R/2 

(47) < | / (i + \yn + i\r 1/2 \G(-y n+1 )\a + k+i|) 1/2 

X J\y„ +1 \>R/2 

x \i[> n \y/n\ipn+i\<iyi ■ ■ -dy n+ i . 
Applying the Schwarz inequality and Lemma f7.lt we obtain from (J4*TI) 

(48) El^l^^CW'+IIWII 9 )- 

n 

Inequalities (|46|) and |48|) complete the proof of Lemma 14.41 □ 
4.1.4. Localization error for Af(0) 2 . 

Lemma 4.5. Let ij) G ^(^(O) 2 ) n D(N f ), and let for some p G (6/5, 2] 

\\V k a x (k)^\\ e L?°(R 3 ) + L 2 (R 3 ). 

Then, the inequality 

(49) (i4 / (0)W,W J ty) + (^(OfV^V'V) - (^/(0)¥,^) < ^ 
ZioZds wzi/i 5 = (j>o — 6/5)/2. 
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Proof. Using the canonical commutation relations, we have 

^(O) 2 = D(0) 2 + D*(0) 2 + 2KeD* (0)L>(0) + cl, 

where the constant c depends on the ultraviolet cutoff. Therefore, it is sufficient to compute 
the localization error for D(0) 2 and D*(0)D(0). We have 

(D(0) 2 U R ip 7 U R ip) + (D(0) 2 V R ip,V R ip) - (D(0) 2 ij,tfj} 

(50) =^Vn + Wn + 2 / G{y n+2 )G{y n+1 )ip n+2 ^n 

n J 

x {u R +2 u R + v R +2 v R - 1) dyi . . . dy n+2 
In the region where max i= i ) ... )n+2 ^ max{|y ri+1 |, |y„ +2 |}, we find 

(51) (u* +2 < + v*_ 2 v* - 1) (!/!,... , y n+2 ) = ((u R ) 2 + (v R ) 2 - 1) (y u ...,y n ) = 0. 
In the region where maxj = i v .. jn+2 \yi\ = \y n +2\ < -R/2, we have 

U n+2(VU Vn+2) = • • • , 2/n) = 1 

and 

V n+2(VU Vn+2) = v R (y 1: ...,y n ) = 0. 

This yields floTj) in that case. Similarly, in the region where max i= i n+2 \yi\ = \y n +i\ < -R/2, 

equation (JoTj) holds. Therefore, in ()50p. it suffices to carry out the integration in the region 
{(y h . . . y n+2 ) | \y n+ i\ > R/2} U {(y u . . . y n+2 ) I \y n +2\ > Let us consider the integral 

in the first region. The other will be treated the same way. We have 



\ / n + 2y/n + l G(y n+2 )G(y n+1 )ip n+2 ij n {u R +2 u R + v R +2 v R - l) dy x . . . dy 



n+2 



(52) -W*j ^y^ 1 + |yn + l|) W |G(y n+2 )|V^TT|Vn| 

xVn + 2|"0„ +2 |(l + |y„+i|) 2,5 dyi • • -dy n+2 . 
Applying the Schwarz inequality and using (J44j) as in Lemma \AA\ we obtain the estimate 

| ^jjn + ly/n + 2 G(y n+2 )G(y n+1 )ip n+2 ip n 

1 



<C 



(u R +2 u R +v R +2 v R -l)d yi . . . dy n+2 
We have 

(D*(0)D(0)U R ip,U R ip) + {D*(0)D(0)V r ^ } V r iIj} - (£>*(0)D(0)V>, V) 

= ^(n + l) / G r (y n+ i)G?(z n+ i)V'n+l(yi > --- > ?/n,yn+l)V'n+l(yi>--- 
n J 

X ( U n+l(Vl> ■ ■ ■ > 2/n> 2/n+l)Wn + l(2/l, • • • , 2/n, Zn+l) 



+ w n+i(?/i' • • • > fn, yn+i)v R (yi, ...,y n , Zn+i) -l)dy 1 ... dy n+ idz n+ i 
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As before, in the region where both y n+ \ and z n+ i are less than R/2, the expression inside 
the integral is zero. Without any loss of generality, we may assume that y n + 1 > R/2. In 
that case, the expression above is bounded by 

(n + l)R 1+s \\iP n+1 ( yi , ...,y n ,y n + l) X (\y n + 1| > R/2)G(-y n + 1)|| 2 

+(n + l)ir ||Vn+i(yi, ...,y n ,z n + l)G(-z n + 1) f 

Similarly to (|52|) . we obtain 

||Vn+l(2/l, • • • ,2/n,2/n + l)x(|y» + 1| > R/2)G(-y n + l)f < fl-^UVrH-lll 2 - 



Therefore, 

(D*(0)D(0)U R ^,U R ^) + (D*(0)D(0)V R ij,V R ij) - (£>*(0)£>(0)V>, V) < j^ s - 
This concludes the proof. □ 
4.1.5. Localization error for the operator a. -8/(0) . 

Lemma 4.6. Let ip G 0(cr.i?/(O)) nD(Nf), and assume that there exists po e (6/5,2], such 
that 

\\V k a x (k)^\\ eL P0 (M 3 ) + L 2 (M 3 ). 

Then, the inequality 

(53) (a.B f {Q)U R i>M R ^) + (o-.B f (0)V R i>,V R i;) - (a.B f {0)^) < 
holds with 5 = (po — 6/5)/2. 

The proof of Lemma 14.61 is similar to the one of Lemma 14.51 with a large number of simpli- 
fications. 

4.2. Proof of Theorem SHI We let 

(54) <S> R :=U R n , 

where Qq is a normalized ground state eigenfunction of the operator T(0), and where IA R is 
defined in (J2Ej). We recall that we have (T(0)Q ,Q ) = £ ||^o|| 2 - We would like to show 
that the value of the quadratic form associated to T(0) at <& R is, for large R, close to the 
value of the quadratic form associated to T(0) at Q . 

First, we notice that Qq fulfills all the conditions of Lemmata 14.1114.61 which implies 

that 

(T(oH,fi > = (T(o)u R n ,u R n ) + (T(o)v R n ,v R n ) + 

ri 

where C{R) tends to zero as R tends to infinity. Thus, since (T{0)V R Q , V R Q ) > So||V R fi || 2 , 
we obtain 

(T(0)$*$*> < E Q + i^l -E (T(0)V^o,V^o) 

< S (1-||V^ || 2 ) + ^ = S ||^|| 2 + M^, 
which proves ii) of Theorem 14.11 
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To complete the proof of Theorem 14.11 it suffices to prove the two Inequalities (J2U|) and 
(1231) . Let us start with (I2U1) . 

\(D(x)^ R ^ R )\<J2^Tl I \G(x-y n+1 )\ \$ R +1 \ \<f> R \d yi dy n+1 

n J 

= + / \G(x - y n+1 )\(l + \x - y n+1 \)\$ R 

\ x \J\y n +i\<R 



R 
+ll 



By applying the Schwarz inequality, we get 

(55) \{D{x)<$> R ,$ R )\ < J^V^TT [\G(x-y n+1 )\ \<$> R +1 \ \® R \d yi dy n+1 

= Enl' G ( x - + \ x ~ yn+i\)(i + \yn + i\r 25 ^ R 



x\\V^TT(l + \y n+ i\) 25 K + u 
We recall that from Lemma 1711 that \G(x - y n +i)(l + \x — y n +i\)\ G L r (M?) for all r > 2. 



Therefore, for p > 3/(3 — 25), and q given by 1/p+l/q = 1, we have || (1 + |?/n+i|) 2S \\q < °o. 



Thus, 

\\G(x - y n+1 )(l + \x- y n+1 \)(l + \y n+1 \y 2S <t> R \\ 

< \\G(x - y n+1 )(l + \x- y n+1 \)x(\ yn+ i\ < R)\\ P \\(l + l3/n+i|) _2 *IUI^||- 

Moreover, for \x\ > 2R, the norm \\G(x — y n+ x)(l + \x — y n+ i\)x(\y n +i\ < R)\\p tends to zero 
as R — > oo. This estimate together with (JH5)) yields 

\(D(x)$ R ,$ R )\ < ^(x)E(||^|| 2 + ||v^TT(l + k + i|) 2 ^ +1 || 2 ) 

Conditions Ciiii) and £\iv ) together with the above inequality conclude the proof of (j2"Uj) if 
we pick S = (po — 6/5)/2. The proofs of (j2"Tj) . (J2~2"j) . and (}2"3*j) are similar. 

5. Approximate ground state for a system with an external potential 

In the present section, we consider the Pauli-Fierz Hamiltonian for M electrons with 
an external potential 

M 

Hm = Yl { H V ^ ® J / + VaA f (x t )) 2 + y/aa ■ B f (xt) + V(x e ) ® J/| 
f=i 

+ 2 ~ x <) ® J / + /e ' ® H / ' 

l<k,e<M 

acting on TC = Hfj <S> T . The brackets ( • , • ) will from here on denote the scalar product 
on 7i. Furthermore, for the rest of this section, we will write operators of the form I e i® Af 
or B e i ® If on TL simply as Aj or B eh respectively, in order not to overburden the notation. 
The precise meaning will be clear from the context. 
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We assume that the Condition £ 2 is fulfilled for this system, which implies, in particular, 
that the operator Hm has a ground state. We will construct an approximation to the ground 
state which is spatially localized with respect to the electron and photon variables, and whose 
energy is close to the ground state energy. 

5.1. Localization of the electrons. We start with localization in the electron configura- 
tion space. To this end, we recall from (EJ) that T denotes the ground state of Hm- For u 
given by ©3, we define T R = {T R , Tf , . . . , T* . . .) e H = H M ® T by 



T n= u \ vT |/ '' |T W , 

where T n is the n-photon component of T. Notice that on the support of T^, we have 
\xi\ < R/2 for % = 1,...,M. 

Lemma 5.1. For all R > 1, 

(56) (H M T R ,T R )<E M + ^ 

(57) i _ < || T *|| < 1 

The proof of this Lemma follows immediately from standard localization error estimates 
for Schrodinger operators [S], and the Condition £ 2 iii). 

5.2. Localization of photons. Our next goal is to localize all photons in a ball of radius 2R 
centered at the origin. For this purpose, we define the function *f? R = (\l/f , ty R , . . . *$f R , . . .) £ 
Hm ® T as 

(58) ^ R = U 2R T R . 

where U R straightforwardly extends the operator defined on T in (|2*Hj) to H M ® 3~ '■ 

We note here that the localization radius for photons is chosen to be four times larger 

than that for the electrons. The consequence is that the contribution of the "external" 

photons to the magnetic vector-potential will be negligible within the region where the 

electrons are localized. 

Similarly to Lemma f4.1[ we find that there exists c < oo, such that for all e > 0, and 

all R large enough, 

(HfU 2R T R , U 2R T R ) + (H f V 2R T R , V 2R T R ) - (H f T R ,T R ) 

Obviously, it suffices to compute the localization error only for the operator 

(~*V X1 + y/aAfix^) 2 + ■ B f { Xl ) + H f . 
In the rest of this section, we will denote x — X\. 
Lemma 5.2. The following estimate holds 

:n \ f ( x )u 2R T R ) + (V 2R T R ,iV x A f (x)V 2R T R ) - (T R ,iV x A f (x)T R )\ 



(60) 



<^(II^ T/? II 2 + Iiv^|| 2 ) 
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where 5 = (p — 6/5)/2 and p is given by £ 2 H)- 

Proof. The proof of this Lemma is very similar to the one of Lemma 14.41 

\(u 2R r R ,iv x D( x )u 2R r R ) + (v 2R r R ,iv x D(x)v 2R r R ) - (r R ,iv x D(x)r R )\ 

(61) < / dxVV^+T f\G x (x-y n+1 )\ \T R +1 \ \V X T R \ 



2 n 

x {u 2R u 2R +l + v^v™! - l)d Vl . . . dy n+1 . 

Similarly to Lemma 14.41 we show that {u 2 ^ i u 2 ^_ 1 + v 2R v 2 n R _ 1 - 1) is nonzero only if \y n+1 \ > R. 
This implies \x — y n +i\ > lz/n+il/2 > R/2. Therefore, the integral in (joT|) can be estimated 



by 

(62) ^l x l^2 n 

x W^Gxix - y n+ i)(l + |x - yn+il) 1 - 5 !! IIV^H 

A 

Since the term ||V X T R || is finite, the rest of the proof is not different from the one of 
Lemma 14.41 □ 

Similarly to Lemmata 14.51 and 14.61 and the above Lemma 15.21 one can prove that 

(63) | (U 2R T R , A\x)U 2R T R ) + (V 2R T R A 2 (x)V 2R T R ) - (T R A 2 (x)T R ) | < 
and 

| (U 2R T R , a ■ B f (x)U 2R T R ) + (V 2R T R , a ■ B f (x)V 2R T R ) - (T R , a ■ B f {x)T R ) \ 



c 

< 



(64) 

* R 1+s 

Theorem 5.1 (Energy of the approximate ground state). For arbitrarily fixed e > and R 
large enough, the following statements hold. 

i) 

(65) E M \\* R \\ 2 <(H M * R ,* R )<E M \\*Y + ^\\* R \\ 2 

R 

ii) Let z G M 3 be an external variable, i.e., the function ty R does not depend on z. Then, for 
\z\ > AR 



(66) \(D(z)* K ,* H )\ < 



(67) \(D(z) 2 * R ,V R )\ < 



c{z) 



\z 



2 



(68) \(D*(z)D(z)S> R ,V R )\< 



c(z) 



\z 



2 



and 

(69) \(K(z)V H ,V H )\< 



where c(z) is a function independent of R that tends to zero as \z\ tends to infinity. 
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Proof. Applying Lemma l5~Tl and Inequalities (|6L)|). (|£>3j) . (J53J) . we obtain 

(70) (H M U 2R T R , U 2R T R ) + (H M V 2R T R V 2R T r ) -L< Em + ^- 
Using EmIIV^T^H 2 < (HjifV 28 ^^^), we get 

(71) (H M U 2R T R ,U 2R T R ) < E M \\U 2R T R \\ 2 + L + 

Since || — > 1 as i? — > oo, we get 

The proof of ii) is analogous to the proof of Lemma f5. 21 and Theorem 14.11 iii). □ 

6. Proof of Theorem 12.11 

The previous discussion enables us to construct a normalized trial function T R,b £ 
H% <g> J". For given JV G N, we define vj/^" 1 as 



R,N-1 



\\&R\ 



where ty R is the function defined in Section 15.21 for a system of M = N — 1 electrons. As a 
natural candidate for a trial state for the proof of Theorem I2.H one could consider the state 
y? = {<p ,(pi, ■ ■ •) defined by 

n 

(72) £j 

x • • • , Vn, Aj+1, . . . , An, X 1} . . . , Xtv_i, Si, ... , sjv-i). 

However, since the components y? n are neither symmetric in the photon, nor antisymmetric 
in the electron variables, our next goal is to symmetrize the function tp n in the photon 
variables, and to antisymmetrize it in the electron variables. 

We denote by S n j the set of (™) possible partitions g of the set of n indices {1, . . . , n} 
into two subsets C\ and C 2 with j and n — j elements respectively. Let i 1 (g), . . . ,ij(g) be 
the indices in C\ and ij + \(g), . . . ,i n (g) in C 2 - We define the function 

(n^Odef Sff^- 1 ) ....//,..a, A n , x, a ) 

\R,b 



' (l/ii, • • • , J/ij, A ii; . . . , Xi,,x N , sjy) 



x ^n-j^iVij+n ■ ■ ■ > 2/«n> A ij+1 , . . . , A in , Xi, . . . , XjV-i, Si, • • • Stf-l) 



Evidently, 



(73) f*» : = y: r) E n u^f 6 *?4 _1 

is symmetric with respect to the permutation of photon variables. 

To construct a combination of the functions T R,b which is antisymmetric in the electron 
variables, let us consider the set of all transpositions 7Tj i — 1, . . . , N, which exchange a pair 
of electron variables (xj, Sj) with (xn,sn), including the trivial transposition (xn,sn) <-> 
(x/v, s^). For an arbitrary function (p(x\, . . . , xn, si, • • • , sjv), let 

(nfV)(xi, . . . ,X N ,Sx, . . . , s N ) := ip(TTi(xi, ...,x N ,s h .. .,S N )). 
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Then, we define 

-1/2 N 



n , x -1/2 JV 

T n, b = ^N-^n J2 E (-i) K(i) nfn^.(,)ef 

.7=0 i=l oesv,.,- 



(74) 

j=o SJ 7 »=1 geSn,. 

where = if i = N, and = 1 otherwise. Obviously, 

r R ' b = (r^ b , rf ' b , . . . , r R > b , ...)en%®T. 

Notice that T R,b is a normalized function in 7i N ® JF, since if |6| > 5i?, the summands in 
(174}) have for different i disjoint supports in electron variables, and thus 

(75) ||r*- 6 || 2 = ||f R ' b || 2 = i 

and 

(H N T R > b } T R ' b ) = (H N t R ' b , f R ' b ). 

Although the state T R,b is not antisymmetric in all electron variables, the quadratic form of 
H N at T R ' b is well-defined. 

Furthermore, both functions Q R,b and vl/^'^ -1 have a finite expectation number of 
photons, say, N\ and N2, respectively. Evidently, this implies that T R ' b has a finite expected 
photon number N\ + A^. 

We remark that for \b\ > 5R, and each of the terms in the sum 

nf E KAtief^tV, 

R ' b and ipn^j' 1 h ave disjoint supports, thus one finds 

(H N T R < b , T R > b ) = (H N f R ' b , f R ' b ), 

where, as we recall from (USD, f R ' b = (f R ' b , f f' b , . . .) is the state prior to antisymmetrization 
in the electron variables. Hence, instead of estimating the quadratic form of the operator 
Hn with respect to the state T R,b , we may estimate it with respect to T R,b . Although this 
state is not antisymmetric in all electron variables, the quadratic form is well-defined. 

We recall that in our notation for the state T R,b , the variables (xn, sn) are the argu- 
ments of Q R > b , while (xi, . . . ,xn-i, Si, • • • , Sjv-i) are the arguments of ■0 ii,iV-1 , and further- 
more, that ||r R ' 6 || = 1. 

Lemma 6.1. For \b\ > 8R, there exists c > independent of R such that the following 
estimate holds 

p3/2 

(H f r R ' b , r R > b ) < (HfQ R ' b , e R ' b ) + (H f y R > N -\ y^- 1 ) + c ^{N f v R \ r R > b ) 

Proof. We have 

n / \ -1/2 

(76) H f f R > b = n\V yi \Y: (" E KM®f' ^nV ■ 

3=0 geS nJ 

Let us start with one of the functions in the sum f)76j) . We take for example the expression 

n| V yi |0 ■ ' (1/1, ... , yj)^f n '_j~ 1 (yj + x, . . . , y n ). All other terms can be treated similarly. 
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In the quadratic form (HfT R,b ,T R,b ), this term appears twice, in 



™(|V yi |©f 6 • • • , yj)^n-j \yj+h ■ ■ ■ , Vn), 
Qf\yi, y j )^^f 1 (y j+lj ...,y n/ / . 



and in 



(77) 



n {\V yi \®j' (3/1,..., (y j+1 ,...,y n ), 
0fl\(2/2, • • • , 1^+1(2/1, Vj+i, -..,!/« 



All other cross terms appearing in the quadratic form (HfT R,b , T R,b ) that contain the function 

n l^2/J©i ' (^1) ■ ■ ■ > %')^n-i _1 (%'+i> • • • > Vn) are zero > because at least for one variable, the 
supports of the functions in the scalar product are disjoint. Let us now estimate (|77jl . The 
function 

©f b (y h • • • , 2/n) 

is supported in the region {\yi\ > |6| — 2i?} whereas 

©f_i(?/2, • • • , Vj+i, ...,y n ) 

is supported in the region {\yi\ < 2i?}. Applying Lemma f7.2l with |6| > 8-R, we arrive at 

Vj©f 6 (1/1, • • .,y j )^^~ 1 (y j+1 , . . .,y n ), 
©f-iG/2, • • • , (j/l, 2/i+i, • • • > 2/n) 

/.. _ D L . ""DAT 1 , 

< C 71: 



(7 



which implies 



| 6 |5/2 (J^ ' ^' ' ' ''V^n-j (%•+!. •••:?/■ 

+ || ©^1(2/2, • • • , %0*^5+f (yi, • • • , 2A. 

^3/2 



2 



(79) 



n / \ -1 

+ E"E (") E (Iv.jnVtojef'er'.nLtojef 1 *^- 1 ). 



n j=0 v 7 se5 n , 

For fixed n and j, in the sum 
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the variable y% appears ( n _J) times in Q R ' b and ( n ™ -^J times in ^f^j 1 . Therefore, the 
second term on the right hand side of (|79|) can be rewritten as 



Ei^)" 1 ^:^^ 



n i=l 



n 

n i=l 



+ EE n (j) (n-j-1 



-i 



)-R,6||2 



),* ft Lj (yi,...y n -j)) 

n 

EE^i^fV • • 8 f • • • >%-)>ii*£? _1 ii 2 



(80) 



n j=l 

n 



+ E E(" - ^')ll©f 1 2 (l V yi l^ 1 ^, . . • yn-j), 9^-\yi, ■ ■ ■ Vn- S )) 

n j=l 

= (H f O R > b , R > b ) + (H f ^ N ~\ V^- 1 ) 

The relations (|75j) and (JHUj) imply the statement of the Lemma. □ 
Lemma 6.2. For any e > and \b\ large enough, 

N 



51) 



]T*v^(*,)f^,f^> 

AT-l 



2{\b\-2R) 

jV-1 



|v xw e^|| 2 + ||e^|| 2 ) 



^2(\b\-2R) V Vxe 11 + " 11 J 



Furthermore, 



N 



(J2vB(x e )t R ' b ,f R ' b ) 
(82) < (a ■ B(x N )Q R > b , <d R > b ) + ^ {° ■ B^)**^" 1 , v^" 1 ) 



2V-1 



=i 



AT-l 



+ — - — iie^ii 2 + V — - — ii^- 1 !! 2 

+ (\b\-2R) 11 11 ^ 2^ (\b\ - 2R) 11 11 
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and 



(83) 



N 



N-l 



(^D 2 (x ^ )f^^f i^ ' f, )-( J D 2 (a; J v)e i^ '^e i^ ' b )-^( J D 2 (x,)^ 7V - 1 ,^^- 1 ) 



N-l 



< 



(\b\-2R) 



ie 



R6112 



E 



? (|6|-2i2) 



I* 



R,iV-l||2 



Moreover, 



N 



1=1 



(J2D*(x e )D(x e )f R \ f R ' b ) - (D*(x N )D(x N )Q R ' b , Q R <") 

N-l 

e=i 

N-l 



< 



(\b\-2R) 



10 



R6112 



2R) 



R,N-1\\2 



Proof. We recall that in . ' - 1 , the variable x\ appears only in the function ^^Lj 1 , 

Ft b 

and the variable x^ only in • ' . Permutations of photon variables do not change this fact. 
We have, for k = 1, . . . , N, 



D{x k )f R > b 

n 



3=0 



n-l 
-1/2 



Vn{G{x k -y n ),Ii p nj (g)Qf b ^ R ^ j l )tfQg?®C*,&Vn) , 



where as before, dy n means integration with respect to y n and summation over the associated 
polarization A n . 

Let us start with one of the functions n£ -(p)Gj^ ^ ' - -1 in the sum (J73J). For fixed 



g, two variants are possible. Either the index n is in C\, and the function Q R,b depends on 
the photon variable y n , or the function ^^f^ -1 depends on y n . For fixed n and j, the first 
variant occurs times, whereas the second one occurs times. Let us consider 

the function 



-1/2 



0j {x N ,y 1 ,...y j )W n _ j (x 1 ,...x N - 1 ,y j+1 ,...y n ) 



r R,N-l 



In the quadratic form (Y^=i iV Xk D(x k )T R ' b , T R ' b ), it appears only once in the scalar product 
with 



\/nV Xk G(x k - y n ) 



n — 1 
3-1 



-1/2 



0f ,6 (x7v,yi, . ..yj) 



XV, i x h ■ ..x N -i,y j+ i, . ..y n -i, 
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which, in the case k ^ N, is equal to 

-l 



(86) 



n — 1 
n-j-1 



(xi, . . .x N -x,y j+ i 



Vn) 



V Xk G(x k - Vn^nLj^ (xi, . . . x N -i,y j+ i, . . . y n -i\ 
and in the case k = N, 

n — 1 



)-R,&||2 



■J 



n-j 



(87) 



n-j- 



x N -i, y j+1 , . . . y n )V XN Qf ,b (x N , yx,... yj) 



G(x N - y n )*M_J(zi, . . . x N - 1 ,y j+1 , . . . y n -i)Qf'°(x N , y u ... y s ) 



R,b, 



All other terms in (J73|) with the same j, and with y n in ^f n '_j , give the same contribution 
to (J2k=i iV Xk D(xk)r R,b , f R,b ). Summing up these ( n Z~li) contributions in the case k ^ N 
yields 



R,N-1 vy 777 7,f,i?,JV-l 



and in the case k = N, 



(89) v '" - /(^V^ef^c^-^^zjef ) 



If we sum first over m = n — j , and then the terms (}88j) over j , we get 
(90) ^ R ' N -\V Xk D(x k )^ N - 1 )\\Q R ' b \\ 2 - 

Let us compute first the sum over n — j of the terms (|89jl. and estimate them according to 
(ffifij) . We obtain for e > 0, and |6| sufficiently large, 



(91) 



R,b\ 



c(x N ) 



R,6||2 



lie 



" |a*| " " 2(|6|-2i2) 

where we used that |xjv| > |&| — 2i2, and c(xn) tends to zero, as \xn\ tends to infinity 
Therefore, 



(92) 



n j 



< 



2(|6| -2/2) 

In analogy to (}90j) and (}92|) . the contribution to (X) fc=1 ^V Xfe D(xfc)f f fl,fe ) of the terms for 
which the variable y n is in Q R,b , is, for k = N, equal to 

(93) (W XN D(x N )G R > b ,e R ' b ) 

and for k ^ N, it can be estimated by 



(94) 2(\b\-2R) 
This completes the proof of (fHTj). 



l|V, fc * 



RiV-l||2 



+ 11* 



R,N-1\\2 
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Let us next prove the inequality (J83j) . The operator D 2 (x k ) acts as 

-1/2 



(95) 



(^) ^ a = E(")" Evw^r 



x - y n )G(x k - y n _i), Yi p n j (g)Q j ' ^(x,^)®!^,^! 



) ' 



where X := M 3 ® C 2 . Assume that in the decomposition g, we have the indices n G C 2 and 
(n — 1) G C 2 . Then, both variables y n and y n _i appear in the function ^fi^ -1 . For fixed 
n and j, we have ( n ™^ 2 ) such cases. Similar to (jHEJ) in the case k ^ N, and to (jHSj) in the 
case k = N, we obtain, respectively, 

,R,N 

' n _j ,^J\^k — yn) ^J\^k — yn-lj ^ n _j_2 



(96) v^jv/n-l-j/^" 1 , G(x fe - y n ) G(x fc - i^O^Zj) ||ef 6 || a 



and 



(97) V" -jy/n-3- l(*2? _1 ©f*, G(a* - y„) G(a* - y^^I^f 

Now, summing each of these expressions over m = n — j and j, and applying ()67|). we arrive 
at 

(98) (D 2 (x k )^ R ' N - 1 ,^ R ' N - 1 )\\Q R ' b \\ 2 
for k N, and 

09) £ lie^ll 2 

for k = N. 

Let us now consider g with n G Ci and (n — 1) G Ci, which implies that the variables 
y n and y n -i are in • ' . We get 

(100) (D 2 (x N )e R '\ © R ' b ) u***- 1 f 

for k = N, and 

(101) ii^- 1 !! 2 

1 J (|6| - 2.R) 2 11 11 

for k^N. 

Finally, let us address the case where one of the indices n, n — 1 belongs to C\ and the 
other one to C%. In this case, one of the variables y n and y n ^\ appears in and the 

other one in 0f\ We have 2 ( n : 2 ) such cases. Note that in each such case, either \G(x k — y n )\ 
or \G(xk — Un-i)\ is small, and the contribution of the sum of these terms can be estimated 

as 



(102) — ((N f Q R '\ Q R ' b ) + (N f * R > N - 1 , V R > N ~ l ) 

|o| — 2R V 



< : ( ||0^ 6 || 2 + \\q,R,N-i\\2 

-\b\-2R [l1 " " 



The estimates fEISjl - flM imply (JSHJ) 

The proof of (|82j) is very similar to the one of (J81|) . □ 
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6.1. Proof. To prove Theorem 12.11 we will show that for suitably chosen parameters R and 
\b\, the trial function T R,b satisfies 

(103) (H N r R >\r R > b ) <e n ^ + e . 



We recall that 

N 



H N =J2{ (-«V XJ + VaA f (x e )) 2 + y/aa ■ B f (x e ) + V(x e )} 
e=i 

(104) +1 W(\x k -x e \) + H f , 

l<k,i<N 

and that, as was shown in the previous section, the inequality ()103jl is equivalent to 

(H N r R > b ,f R > b ) <E N ^ + E . 

For M G N, we define 

m 1 

I M (x 1 ,...,x M ) = ^2v(x e ) + - ^2 W(x k -xt). 

1=1 l<k,£<M 

Obviously, we have 

JV 

X)<-A«f R \ r R > b ) + (i N ( Xl , x n )f R > b , r R > b ) 
t=\ 

JV-1 

(105) tr 

+ (i N ^(x u . . . , a*-!)**-"- 1 , ^ N - 1 ) + (-Ae R >\ e R > b ) 

/ JV-1 

+ ( v(x N ) + (J2 w(x t - x N )* R ' N -\ *^-i>e^, e* 6 
\ i=i 

where we used that Q R,b and xl^^ -1 are normalized. On the support of the function Q R ' b , 
we have \xn\ < + R and on the support of the function xj/^^ -1 , \xi — a;Ar| > |6| — 2R. 
This implies, for sufficiently large, that on the support of T R,b (defined in (HHJ)), 

(106) %) + i;Wfe-^)<-^ + 2 i^< " 



1=1 



|6| +i? |6| -2i2 2|6| 



for v = 70 - 7i (iV - 1) > 0. Thus, (fTUHj) and (fTO yield 

^(_A,f^, f + (^(xx, . . . , x n )f^ fe , f^ fe ) 

£=1 

JV-1 

(107) < ^(-A,^- 1 , ^f^- 1 ) + (/jv^l, . . . , X N ^ R ' N -\ V^- 1 ) 



+ (-Ae R > b ,e R > b )-^. 
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Taking into account that HV^.^ 



R,N-1\ 



< c||W 



R,N-1\ 



(£ = 1,. . .N — 1), and that 



I V Xjv R ' b || < c||© ' ||, with a constant c independent of R, we derive from ()81|) 



(108) 



A/ 



Ar-1 



£=1 



(V XN A(x N )e R > b ,Q R > b ) 



< 



\b\-2R 



Similarly to l|I0g|) . and using (J2U, (jH7 |l . and (fBBjl . we have 



A' 



^(A 2 (x,)f^,f^) 



(109) 



£=1 



\b\-2R 



Along the same lines, we have for the magnetic term, using (f2*H|) and (jEHJ), 



N 



(110) 



N-l 



< ■ B(x e )* R ' N -\ ^ R ' N - 1 ) + (a • B(x N )Q R < b , e R > b ) + —3 



|6| -2R' 



According to Lemma f6. II we have 

(H f f R ' b ,f R ' b ) < (H f e R ' b ,e R ' b ) + (Hf***"- 1 ,**"- 1 ) 



Equality ((75J) implies that (N f T R < b ,T R < b ) < c (jl^^- 1 !) 2 + \\Q R > b \* 2 

Collecting the estimates f)107|) - f)l 1 lj) we obtain for any e > and sufficiently large R, 



(H N T R > b , f R > b ) < e n ^ + s - — + 



6c 



+ 



zR?l 2 



2\b\ \b\-2R |6| 5 /2 ' 

To complete the proof of the Theorem, we pick first R large enough to have e < 48~V, and 
then pick |6| sufficiently large to satisfy the inequality (i?|6| -1 ) 3 / 2 < 5(4c) _1 , which implies 

(H N r R > b ,r R > b ) <e n ^ + e . 



Lemma 7.1. We define G\ as 



7. Appendix 



where T denotes the Fourier transform. Then, for A = 1, 2 and arbitrary e > 0, \G\(y)(l + 

M)|eL 2 +^(M 3 ). 
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Proof. The statement of the Lemma follows from the Hausdorff- Young inequality, and the 



fact that for arbitrarily e > 0, 



is in L? e (M 3 ), for i = 1,2,3, which can be 



checked directly. □ 

Lemma 7.2. Let <f\{x) E iJ 1 / 2 (IR 3 ) with support in the ball of radius aR centered at the 
origin, and <p 2 (x) E ff 1 / 2 (lR 3 ) with support outside the ball of radius bR centered at the 
origin. Then for b > 2a, 

1 fl 3 / 2 

(H2) Kiv^.^i^^^-^a^f+feR 

Proof. Consider the function u defined in (12*7)) . Then, for Xi( x ) = u {\ x \/(bR)) an d X^{ x ) — 
\J\ — we have, according to [TUl Theorem9] 

(|V|((^l + <p 2 ),<Pl + V2) - <| V|V?1, V?i> - (|V|^ 2> V2) 



< 



2tt 2 




\ipi(x) + <p 2 (x)\ 


<Pi(y) + <P2(y)\ 




x — 


y 


4 



Yl \xK x ) -xKy)\dydy 



i=l,2 

Since Xi — 1 011 the support of cpi, Xi — on the support of <f 2 , we obtain 

<|V|(v?i + ^2), Vi + ^2) - (|V|v3i, y>i) - (|V|v? 2 , P2) 
= 21Ze{\V\(p 1 ,<p 2 ) 

!_ I I I^WI 1^)1 

^ J J f - y| 4 

/ll 1 1 2 11 n2\ 

" 7r3V2i2(6- a )5/2( M +l|(/9211 } 



□ 
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